Multifractality of wavefunctions at the quantum Hall transition revisited 
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We investigate numerically the statistics of wavefunction amplitudes ip(r) at the integer quantum 
Hall transition. It is demonstrated that in the limit of a large system size the distribution function 
of \ip\ 2 is log-normal, so that the multifractal spectrum f(a) is exactly parabolic. Our findings lend 
strong support to a recent conjecture for a critical theory of the quantum Hall transition. 
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In 1980 von Klitzing, Dorda, and Pepper discovered 
that the Hall conductance cr xy of a two-dimensional elec- 
tron gas develops plateaus at values quantized in units 
of e 2 /h. Although two decades have passed, the integer 
quantum Hall effect (IQHE) still constitutes one of the 
great challenges of condensed matter physics. Initially, 
the effort focussed on the physics of the Hall plateaus 
which is by now fairly well understood Then in- 

terest has shifted towards the transition region, where 
cr xy crosses over from one plateau to the next. However, 
here the situation is not as well resolved. It has been un- 
derstood from early on that this is a second-order phase 
transition, and the scaling scenario has been confirmed in 
numerous experiments and computer simulations [Q. In 
particular, the numerical prediction for the critical index 
of the localization length v ~ 2.35 has been verified ex- 
perimentally. In contrast, analytical approaches did not 
lead to quantitative predictions. The ultimate goal here 
is to identify the effective low-energy theory of the crit- 
ical point (expected to be a conformal field theory) and 
to calculate, on this basis, critical exponents and other 
characteristics of the transition region. 

The earliest field-theoretical formulation of the prob- 
lem was given by Pruisken Q (see Q for a more precise, 
supersymmctric version) and has the form of the non- 
linear cr-model with a topological term. Since the latter is 
invisible in perturbation theory, one has to resort to non- 
perturbative means in order to address the critical behav- 
ior. Pruisken and coworkers were thus led to the dilute 
instanton gas approximation 0. However, this approxi- 
mation can only be justified in the weak-coupling limit, 
&xx/(e 2 /h) 1, and becomes uncontrolled in the critical 
region a xx ~ e 2 /h. For this reason, no quantitative pre- 
dictions for critical properties have been made within this 
approach. Another line of efforts was based on a map- 
ping of the low-energy sector of Pruisken's model onto 
an antiferromagnetic superspin chain The superspin 
chain was also obtained by starting from the Chalker- 
Coddington network model of the IQHE Q. However, 
attempts to find an analytical solution of the superspin 
chain problem remained unsuccessful. 

Recently, two papers appeared which may signify a 



breakthrough in the quest for the conformal critical the- 
ory of the IQHE. Zirnbauer |ll]] and, a few months later, 
Bhaseen et al. |ll| , proposed that this is a cr-model with 
the Wess-Zumino-Novikov-Witten (WZNW) term T, 

s[g] = J d2x strcV&r 1 ^ + fcr, (1) 

where g belongs to a certain supersymmetric superspace 
(see [ |10|Jll|1 for a detailed exposition), and A and k are 
coupling constants. The theories considered in [[To[ and 
|TTf differ only in the value of the constant k in front 
of the WZNW term, k = 1 in ffij and k = 1/A 2 in 
pl| . If true, the conjecture of Refs. |lO|jll[| will provide 
a framework for a systematic study of the IQHE critical 
behavior. While the exponent v has not yet been found 
on this basis, a non-trivial prediction for the statistics 
of critical eigenfunctions has been obtained. Specifically, 
it was found in (jnj that the corresponding multifractal- 
ity spectrum is exactly parabolic (which means that the 
distribution of eigenfunction intensities is log-normal), 

f(a) = 2 - (a - a ) 2 /4(«o - 2), (2) 

with ctQ — 2 = 2A 2 (a formal definition of the function 
f(a) is given below). For the model of [[[o| a parabolic 
multifractal spectrum of a somewhat different quantity, 
the two-point conductance, was obtained JT3] . We can 
show, however |l3) that the two statements are closely 
related, so that the findings of Jl(]|l^] imply again the re- 
sult (||) for the eigenfunction statistics, with a — 2 = 4A 2 . 
The exact parabolicity (||) of the multifractality spectrum 
constitutes a stringent prediction, numerical verification 
of which would be a serious check of validity of the theory 
proposed in ||l0| , |Tl"|| . This is the aim of the present paper. 

Our interest in this problem was additionally moti- 
vated by the fact that the previously published numer- 
ical results P,^,^4]-p7t appeared to be in clear conflict 
with the prediction (0), showing strong deviations from 
parabolicity. Furthermore, it has been a widespread be- 
lief that the parabolic approximation cannot in principle 
be exact, since the singularity spectrum f(a) is only de- 
fined on an interval a_ < a < aq_, where it is positive, 
and has infinite derivatives at the termination points, 
/'(<*_) = f'(a + ) = oo. As we demonstrate below, these 
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earlier findings were crucially affected by finite-size ef- 
fects and by the absence of a proper ensemble averaging. 

To begin with, we recall the general framework used 
for the description of critical wavefunction statistics. One 
introduces a set of inverse participation ratios |Ts{ ] 



P q = 



d d r\iP{r)\ 2q , q>0, 



(3) 



where i/j(r) is the wavefunction amplitude and L the sys- 
tem size. (For the IQHE dimensionality is d — 2.) These 
moments show in the large L-limit scaling behavior char- 
acterized by a set of exponents, which can be defined for 
both the average and the typical value of P q , 



{P q )=H q L- f «, 

exp(lnP 9 ) 



ptyp 
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c q L~ 



(4) 
(5) 



where (...) denotes ensemble averaging. As long as q is 
small enough, q < q c , the distribution function of P q is 
sufficiently narrow and f q = r q . However, with increas- 
ing q this distribution function becomes broader. Most 
importantly, it shows a power-law tail oc P q 9 at large 
Pq, where x q decreases with increasing q. The critical 
value q c is determined by the condition x qc = 1. For 
q > q c the average (P q ) is governed by rare events, and 
r q > f q . We refer the reader to Ref. jl5| for more details. 

Instead of using the moments P q one can study directly 
the distribution V of the wavefunction amplitudes. It is 
easy to see that translates, in the limit L — > oo, into 

V(a) ~ (lnL) 1 / a L-* f ^, a = - In \ij)\ 2 / InL, (6) 

where the function /(a) is related to the exponents f q 
via the Legendre transformation, 



fg = qa-f(aj; q = f'(a). 



(J) 



Similarly, one can define the /(a)-spectrum of a typical 
eigenfunction (which satisfies /(a) >0 §), 



f(a) = qa -T q ; q = f'(a). 



(8) 



The two definitions coincide, /(a) = /(a) > [Esl, for 
q < q c , or, equivalently, for a_ < a < ao, where ao 
describes the scaling of the typical value, exp(ln \4>\ 2 ) °c 
L~ a ° , and a_ denotes the zero of/, /(«-) = fipc-) = 0. 
On the other hand, in the regime of rare events, a < a_, 
the function /(a) becomes negative, while /(a) is not de- 
fined. The best estimate for ao from previous numerical 
work is 2.26 ± 0.01 JTeJ . The value of a_ obtained in Jl|] 
is 1.1 ± 0.1, incompatible with (j^). 

It is important to realize that the field-theoretical pre- 
diction (|^) refers to the function /(a), since in the theory 
one deals with averaged moments (^). In contrast, the 
earlier numerical studies were devoted to the multiffac- 
tal spectrum of a single eigenfunction, thus yielding f(a) 



as an output. We further turn to a common misconcep- 
tion concerning an infinite slope of the f(a) spectrum 
at its termination point, /'(a_) = oo, which seems to 
rule out the parabolic form (|^) . Indeed, in any finite sys- 
tem r q must be an analytic function of q, implying that 
/(a) meets the a-axis with infinite slope This is, 
however, nothing but an artefact of a finite system size 
L. For any given L there is a finite interval of a close 
to a_ where the corresponding spectrum /z,(a) deviates 
from its asymptotic form and approaches the axis with 
an infinite derivative. In the limit L — > oo this interval 
shrinks to zero, and the spectrum acquires its limiting 
form /00(a) = /(a), corresponding to a non-analytic T q , 



'q, 

qa. 



0<q<q c 
q > q c - 



(9) 



The situation is similar to a phase transition where sin- 
gularities also occur in the thermodynamic limit only. 
Practically speaking, this means that r q and f q , though 
identical in the thermodynamic limit on the interval 
< q < q c , have different finite-size corrections. 

For calculating the critical wavefunctions at the quan- 
tum Hall transition we employ the Chalker-Coddington 
network model f2(i|j . In order to obtain the wavefunction 
we translate the lattice dynamics into a unitary time evo- 
lution operator U which describes the wave packet prop- 
agation on the network in discrete time steps |fl6|| . The 
desired critical wavefunctions are the eigenfunctions of 
U, which are found by numerical diagonalization. 

The calculations were performed for systems of the size 
ranging from L = 16 (6 • 10 6 wave functions) to L = 1280 
(2000 wave functions). We implemented efficient numeri- 
cal packages |ji"|-p3| which allowed us not only to (partly) 
diagonalize large systems, but also to do it fast enough 
in order to collect sufficient statistics. Specifically, for a 
system of size L and for each disorder realization we di- 
agonalized a complex N x N matrix of the size N = 2L 2 
(which reached N w 3.3 • 10 6 for the largest L) with ma- 
chine accuracy. For comparison, the biggest realization 
of the network model reported in the literature that we 
are aware of is smaller by a factor of five in linear di- 
mensions Jl6| . Since an accurate extrapolation to the 
thermodynamic limit was of primary importance for our 
work, it was crucial that we could observe the finite size 
corrections over almost two orders of magnitude in L. 

As a first test of the validity of Eq. (g), we check its 
two particular implications, namely ao — 2 = 2 — ai and 
a i/2 = 2, where a q is determined by (0), a q — df q /dq. 
We evaluated a q (L) for different system sizes, 



a q (L) = -m 2 nn\rP\ 2 )/(\^)lnL, 



(10) 



and then extrapolated to the infinite system size, a q — 
lim£_>oo otq(L), according to 



(ii) 
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The most important finite-size effect is of the form 1/ In L 
according to (Q), with a coefficient \ = d\nc q /dq. The 



second term in brackets in (|llj) stems from finite-size cor- 
rections to and is governed by the leading irrelevant 
scaling exponent y. The corresponding length scale £; rr 
can be in principle absorbed in the coefficient j q . 
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FIG. 1. Finite size effects for the scaling exponents ao (°, 
fit solid line) and Qi (o, fit dashed line) defined in (|lo|). Also 
shown is the result for qq (□ , fit dotted line) from the transfer 
matrix calculation, Eq. fll2[). 

The raw data that provide the basis for the extrap- 
olation to L — > oo are shown in Fig. [l| From the fit 
([fl]) we find a = 2.261 ± 0.003, A = 0.09 ± 0.02 and 
ai = 1.739 ±0.002, Ai = -0.11 ±0.01. A similar analysis 
for the case 9 = 1/2 yields a 1/2 - 2 = 0.0001 ± 0.0002, 
X 1 / 2 = -0.042 ± 0.002. It is clearly seen that within the 
accuracy of the numerical data ckq — 2 and 2 — a\ indeed 
coincide and ct\/2 — 2 vanishes, in agreement with (|^). 
The irrelevant exponent y is found to be y — 0.4 ±0.1, 
which is consistent with earlier results (the accu- 

racy of its determination is not very high, since the cor- 
responding finite-size correction is rather small). 

As an additional check we have calculated ao in yet 
another way, which uses a conformal mapping to a quasi- 
one-dimensional strip of the Chalker-Coddington net- 
work. It was shown by Janfien that ao is related 
to the localization length £j, in a strip of width L: 



a -2 = (1/tt) lim L/£ L . 

L — >oo 



(12) 



Since this relation ties ao to the amplitude A of the scal- 
ing law £jr = LA, the finite-size corrections on the strip 
take a form different from (fiOj), 



a (L)=a + j(&„/L)y + 



(13) 



We have calculated £x, using the transfer matrix method 
and found ao — 2.260 ±0.003, in full agreement with the 
previous result. This confirms the fundamental assump- 
tion that the critical theory is conformally invariant. The 



irrelevant scale index is obtained as y 
agreement with earlier findings. 



0.45 ± 0.1, in 




FIG. 2. Evolution of f(a), as calculated from the distribu- 
tion function (|^), with the system size L. The inset shows 
the corresponding change of a_, defined by /(«-) = 0. An 
extrapolation to infinite L yields: a_ = 0.8 ± 0.02. 



Having satisfied ourselves that numerical results for 
ao, oti/2, and a\ are in favor of the conjecture (||), we 
turn to an analysis of the /(a) spectrum in a broader 
range of a. Results of the calculation of f(a) according 
to Eq. (H) from the distribution function V(\ip 2 \) for dif- 
ferent system sizes are shown in Fig. 0. It is seen that 
f(a) changes appreciably with L, which points to the 
importance of the extrapolation to L — > oo. In partic- 
ular, the inset of Fig. || shows the evolution of the zero 
a_ of /(a) with L. Extrapolating to infinite system size 
(by taking into account the leading finite-size correction, 
which has the same 1/lnL form as in Eq. (pT|)) yields the 
value a_ = 0.8 ± 0.02 in the thermodynamic limit. This 
agrees with the result 0.816 ± 0.004 based on Eq. (j|) in 
combination with the value of ao found above, providing 
a further strong support to the parabolic law (0). 

In order to obtain the whole /(a) curve in the thermo- 
dynamic limit with highest possible accuracy, we return 
to the procedure based on the evaluation of moments, 
Eqs. (^), (jl(i|), since it has the advantage that we have full 
control over finite-size corrections, Eq. (11) (the trans- 
formation from the moments to the distribution func- 
tion Eq. (^|) induces additional corrections proportional 
to higher powers of 1/ In L) . Performing an extrapola- 
tion to L — > oo, we get the /(a) curve shown in Fig. [jj, 
which represents the central result of this Letter. The 
figure demonstrates a perfect agreement between the ob- 
tained /(a) and the parabolic form (||) with a = 2.262. 
Of course, a numerical analysis can never guarantee that 
two quantities are identical rather than just very close 
to each other. However, taking into account that there 
is no small parameter in the problem, we believe that 
an accidental closeness of f(a) to a parabola with such 
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a high accuracy is extremely improbable. We thus con- 
clude that the theoretical prediction (||) stating that the 
multifractal spectrum is exactly parabolic is correct. 
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FIG. 3. Multifractal spectrum f(a) (solid lineL The 
dashed line shows the theoretical conjecture (m) with 
ao = 2.262. The point where /(a) vanishes according to the 
data presented in Fig. Q is also shown (□). Inset: spectrum 
/(a) of a typical eigenfunction for different L and extrapo- 
lated to L —* oo (•). The line shows the parabolic law. 

To illustrate the difference and the relation between 
/(a) and /(a) and to make a closer contact to earlier 
works, we show in the inset of Fig. ||the f(a) spectrum 
as obtained for different system sizes. The extrapolation 
to L — > oo gives a good agreement with the parabola 
at a > ao, except for a close vicinity of a_, where the 
finite-size corrections to f(a) become especially strong. 

Finally, we mention an open issue concerning the val- 
ues of the coupling constant A and k in Eq. (Q) deter- 
mining ao in the theory. While Zirnbauer JTy] gave con- 
vincing arguments in favor of A: = 1, the mechanism of 
fixation of A remains unclear. 

In conclusion, we have studied the statistics of crit- 
ical wave functions at the IQHE transition point. We 
implemented a powerful algorithm which allowed us to 
reach unprecedented large system sizes and to gather suf- 
ficiently good statistics. Having performed the ensemble 
averaging and an analysis of the finite-size corrections, 
we calculated the multifractality spectrum f(a) of crit- 
ical wavefunctions in the thermodynamic limit L — > oo 
for 0.5 < a < 2.5. The result is in perfect agreement 
with the parabolic form (||) with a$ = 2.262 ± 0.003 and 
hence it supports the recent conjecture for the form of the 
conformal field theory of the IQHE critical point |l0|,[ll| . 

We are thankful to X.S. Li and to the Computer 
Centers of the Forschungszentrum and the University 
of Karlsruhe for technical assistance. Discussions with 
M. Janfien, R. Klesse, D.G. Polyakov, L. Schweitzer, 
A.M. Tsvelik, P. Wolfie and M.R. Zirnbauer are grate- 
fully acknowledged. This work was supported by the 



[4] 
[5] 
[6] 
[7] 
[8] 
[9] 



[10] 

[11] 
[12] 

[13] 



[14] 

[15] 
[16] 

[17] 
[18] 



[19] 

[20] 
[21] 

[22] 
[23] 
[24] 



Also at Petersburg Nuclear Physics Institute, 188350 
St. Petersburg, Russia. 

K. von Klitzing et al., Phys. Rev. Lett. 45, 494 (1980). 
R.E. Prange, S. M. Girvin (eds.), The Quantum Hall Ef- 
fect (Springer, New York, 1992). 

M. Janfien, O. Viehweger, U. Fastenrath, and J. Hajdu, 
Introduction to the Theory of the Integer Quantum Hall 
Effect (VCH, Weinheim, 1994). 
B. Huckestein, Rev. Mod. Phys. 67, 357 (1995). 
A.M.M. Pruisken, Nucl. Phys. B 235, 277 (1984). 
H.A. Weidenmiiller, Nucl. Phys. B 290, 87 (1987). 
H. Levine et al., Nucl. Phys. B 240, 30, 49, 71 (1984). 
M.R. Zirnbauer, Ann. Physik 3, 513 (1994). 
D.H. Lee, Phys. Rev. B 50, 10788 (1994); J. Kondev and 
J.B. Marston, Nucl. Phys. B 497, 639 (1997); M.R. Zirn- 
bauer, J. Mat h. Phys. 38, 2007 (1997). 



M. Zirnbauer 
M.J. Bhaseen et al 
M. Janfien et al. 



hep-th/9905054| v2 



Nucl. Phys. B 580, 688 (2000). 
Phys. Rev. B 59, 15836 (1999). 
The connection is established as follows. Using the re sults 
of R. Klesse and M. Zirnbauer, cond-mat/0010005, it is 
easy to show that (T/J ~ L d {\f m \' 2{1 - q >\i() l \' 1 ' 1 ), where I 
and m are two points of the network, T; m is the corre- 
sponding conductance, ipi and ip m the wavefunction am- 
plitudes, and q < 1/2. This leads to the following relation 
between the exponents X q characterizing the scaling of 
moments of the two-point conductance, (T^ m ) ~ L~ Xq , 



X a 



(Tim/ 
d + Tq + fl- 



The 



and the indices f„ defined in 

parabolic form (fe|) of the wavefunction multifractality 
leads thus to a parabolic spectrum for the conductance 
too ? X a — Xtq(l — q) (which is exactly what was found 
UTjpflliJ), with X t = 2(q - d). The result X t = 8A 2 of 
[110 1 translates in this context into ao — 2 = 4A 2 . 
B. Huckestein, B. Kramer, L. Schweitzer, Surf. Sci. 263, 
125 (1992). 

I. Varga et al, Europhys. Lett. 36, 437 (1996); 

R. Klesse and M. Metzler, Europhys. Lett. 32, 229 

(1995); Int. J. Mod. Phys. C 10, 577 (1999) 

M. Janfien, Int. J. Mod. Phys. B 8, 943 (1994). 

We restrict ourselves to q > (corresponding to a < ao) 

in this paper. One can generalize the whole consideration 

to the range of negative q (a > ao), but this requires 

some modification of the definition fel). 

F. Evers and A.D. Mirlin, Phys. Rev. Lett. 84, 3690 

(2000); A.D. Mirlin and F. Evers, Phys. Rev. B 62, 7920 

(2000). 

J.T. Chalker and PD. Coddington, J. Phys. C 21, 2665 
(1988). 

PR. Amestoy et al. Comput. Methods in Appl. Mech. 
Eng. 184, 501 (2000); PR. Amestoy et al, SIAM J. Ma- 
trix Anal. Appl. 23, 15 (2001). 

J. W. Demmel et al., SIAM J. Matrix Anal. Appl. 20, 
720 (1999). 

R.B. Lehouca, D. Sorensen and C. Yang, ARPACK Users 
Guide (SIAM, Philadelphia, 1998). 

F. Evers and W. Brenig, Phys. Rev. B 57, 1805 (1998). 



4 



